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Abstract 

A model is proposed for the origin of cosmic rays (CRs) from ~ 10 14 eV/nucleon to 
the highest energies 10 20 eV). GRBs are assumed to inject CR protons and ions 
into the interstellar medium of star-forming galaxies — including the Milky Way — 
with a power-law spectrum extending to a maximum energy ~ 10 20 eV. The CR 
spectrum near the knee is fit with CRs trapped in the Galactic halo that were accel- 
erated and injected by an earlier Galactic GRB. These CRs diffuse in the disk and 
halo of the Galaxy due to gyroresonant pitch-angle scattering with MHD turbulence 
in the Galaxy's magnetic field. The preliminary (2001) KASCADE data through the 
knee of the CR spectrum are fit by a model with energy-dependent propagation of 
CR ions from a single Galactic GRB. Ultra-high energy CRs (UHECRs), with ener- 
gies above the ankle energy at ^ 3 x 10 18 eV, are assumed to propagate rectilinearly 
with their spectrum modified by photo-pion, photo-pair, and expansion losses. We 
fit the measured UHECR spectrum assuming comoving luminosity densities of GRB 
sources consistent with possible star formation rate histories of the universe. 

For power-law CR proton injection spectra with injection number index p <; 2 
and low and high-energy cutoffs, normalization to the local time- and space-averaged 
GRB luminosity density implies that if this model is correct, the nonthermal content 
in GRB blast waves is hadronically dominated by a factor ~ 60-200, limited in 
its upper value by energetic and spectral considerations. Calculations show that 
100 TeV - 100 PeV neutrinos could be detected several times per year from all 
GRBs with kilometer-scale neutrino detectors such as IceCube, for GRB blast-wave 
Doppler factors 5 ^ 200. GLAST measurements of 7-ray components and cutoffs 
will constrain the product of the nonthermal baryon loading and radiative efficiency, 
limit the Doppler factor, and test this scenario. 
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1 Introduction 



In this paper we develop a model for high energy cosmic rays (HECRs; here 
defined as ^ 10 14 eV/nucleon CRs), based on the underlying assumption that 
CRs are accelerated in the relativistic and nonrelativistic shocks found in 
GRBs and their attendant supernovae (SNe). This model extends the work 
of Vietri (0) and Waxman (2), who proposed that UHECRs originate from 
GRBs (see also Ref. (j3jVh to include cosmic ray production from SNe and GRB 
sources in our Galaxy Q). To test the model, KASCADE HH0) and HiRes-I 
and HiRes-II Monocular data (0) are fit over the energy range ~ 2 x 10 16 eV to 
w 3 x 10 20 eV. A good fit to the entire data set is possible if CRs are injected 
with a power-law spectrum with number index p = 2.2, as could be expected 
in a scenario where particles are accelerated by relativistic shocks f)9t llOl). 



CRs that are injected with energies 10 19 eV diffuse through and escape from 
their host galaxy The sources of high-energy CRs, namely GRBs, are located 
in star-forming regions found in the galaxy's disk. CR transport in the disk 
and halo of the Milky Way is modeled using a time-dependent, spherically- 
symmetric propagation model that employs an energy-dependent, spatially- 
independent diffusion coefficient. The random- walk pathlengths are assumed 
to arise from gyroresonant, pitch-angle scattering of CRs with a magnetohy- 
drodynamic (MHD) turbulence spectra that can be decomposed into two com- 
ponents reflecting different power-law distributions of turbulence over different 
wavelength ranges. The spectral break at the CR knee energy is explained in 
an impulsive, single-source model for HECRs if the turbulence spectral index 
changes from Kraichnan to Kolmogorov turbulence near the wavenumber res- 
onant with CRs at the knee of the CR spectrum. We model the preliminary 
KASCADE data reported in 2001 by assuming that all ionic species have 
the same injection index, with the compositions of the ions adjusted to fit 
the observed CR spectrum near the knee. Additional hardenings of the low- 
energy CR spectrum from a single GRB source of HECRs can result both from 
energy-dependent diffusion and from a low-energy cutoff in the CR injection 
spectrum. Superposition of the contributions from many SNe are assumed to 
accelerate the bulk of the GeV/nuc - TeV/nuc CRs, as in the conventional 



scenario ( lit |1 



UHECRs have such large gyroradii and diffusion mean free paths that they are 
assumed to escape directly from the halo of the GRB host galaxy and stream 
into metagalactic space. The UHECR energy spectrum evolves in response to 
photo-pair and photo-pion losses on the redshift-dependent cosmic microwave 
background radiation (CMBR), which we treat in a continuous energy- loss 
formalism. UHECRs in intercluster space also lose energy adiabatically due 
to cosmic expansion. The measured UHECR spectrum arises from the con- 
tributions of sources throughout the universe, with an intensity that depends 
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on the local luminosity density of GRB sources and the evolution of the GRB 
luminosity density with redshift (see Refs. (H3) and (Q) for recent reviews 
of UHECRs and GRBs, respectively). The UHECR spectral model assumes 
that many sources produce the measured ultra-high energy and super-GZK 
(> 10 20 eV) CRs. 

Our model for UHECRs from GRBs implies a local time- and space-averaged 
CR luminosity density Eqr = fcR&GRB,xh °f CRs. The local (z C 1) luminos- 
ity density Egrb^x/^ ~ 10 44 ergs Mpc~ 3 yr _1 is inferred from BATSE observa- 
tions of the hard X-ray/soft 7-ray (X/7) emission from GRBs (flit llil 13). The 
value of Eqr depends sensitively on the minimum energy E min of CR injection 
for soft injection spectral indices p ^ 2. For p ~ 2.2, if E min = 1 GeV then 
^ 700 x more energy must be injected in nonthermal hadrons than is observed 
as X/7 emission from GRBs. If E m i n « 100 TeV, then « 70 x more energy is 
required. Such large baryon loads would provide a bright cascade emission sig- 
nature in GRB spectra at MeV - GeV energies through photopion or hadron 



synchrotron processes (|17t |1 



The large nonthermal baryon load, fcR 3> 1, required to fit the UHECR 
data assuming that the GRBs comoving luminosity density traces the star 
formation rate (SFR) history, implies that GRBs can be much more luminous 
neutrino sources than predicted under the standard assumption that the en- 
ergy injected into a GRB blast wave in the form of CRs is equal to the energy 
inferred from X/7 fluence measurements of GRBs (flih. We therefore predict 
that if this model is correct, and CRs are accelerated and injected in the form 
of soft power- laws with index p ^ 2, then IceCube should detect up to several 
neutrinos from the brightest GRBs with total X/7 radiation fluence at levels 
3> 10~ 4 erg cm -2 . This prediction holds both in a collapsar scenario when 
the Lorentz factors V of the relativistic outflows are ^ 200, or if the GRB 
takes place in an intense external radiation background for a wide range of V. 
Lower limits to V can be inferred from 7-ray transparency arguments applied 
to observations of GRBs with the Gamma-ray Large Area Space Telescope 1 . 



We also consider whether this model can explain the AGASA data (J20J) for 
the UHECRs. Poor fits are found if GRBs inject soft CR spectra with p ^ 2. 
However, if GRBs inject hard spectra with p w 1, for example, through a 
second-order relativistic shock-Fermi process (j2lh or through the converter 
mechanism (22), then the highest-energy AGASA data can be fit, though the 
reduced \ 2 °f our best fits are not compelling. Because the injection spectrum 
is so hard, most of the produced high-energy neutrinos are too energetic and 
the flux too weak to be detected with IceCube, though other telescope arrays, 
such as the Extreme Universe Space Observatory (EUSO), could be sensitive 
to these GZK neutrinos. To explain the CRs below ~ 3 x 10 19 eV, an additional 
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component of CRs would still be required either from GRBs or another class 
of sources, and these would make an additional contribution to high-energy 
neutrino production. 

Section 2 gives a discussion of the CR and photon luminosity density of GRBs 
and their event rate. Our propagation model describing CR diffusion in the 
disk and halo of the Galaxy is presented in Section 3, where we fit the KAS- 
CADE data between w 0.8 and 200 PeV with a single Galactic GRB source 
~ 500 pc away that took place around 200,000 yrs ago. In Section 4, we de- 
scribe our calculation of the UHECR flux, including energy losses from cosmic 
expansion, and photo-pair and photo-pion production. We present minimum 
Xr fits to the high-energy KASCADE, HiRes-I and HiRes-II Monocular data 
covering the energy range 2 x 10 16 eV to 3 x 10 20 eV. We also fit the AGASA 
data from 3 x 10 19 eV to 3 x 10 20 eV for hard CR injection spectra. Section 5 
presents new high-energy neutrino calculations from hadronically dominated 
GRBs, and our predictions for km-scale high-energy neutrino telescopes such 
as IceCube or a deep underwater, northern hemisphere array. Discussion of the 
results and conclusions are given in Section 6. Our treatment of the UHECR 
attenuation and flux calculation is described in Appendix A. 



2 Gamma Ray Burst Model for High-Energy Cosmic Rays 



The progenitor sources of HECRs are likely to be GRBs, as suggested by the 
following observations: 



(1) The evidence from KASCADE (7) that the break energies of the different 
CR ionic species are proportional to rigidity, and that the mean atomic 
mass increases with E through the knee region, admits a propagation 
solution to understand the energies of the spectral breaks of different CR 
ions — given a galactic source that injects power-law CRs to the highest 
energies. 

(2) Composition changes above the "second knee" at E m 10 17,6 eV, and 



possible GZK features in the UHECR spectrum (|23t |24J), implies that 
UHECRs and their sources are metagalactic. 
(3) The CR all-particle spectrum breaks at = 3 PeV by « 0.3 units and 
then extends without spectral change, other than for a possible weak 
softening above the second knee, to the ankle at E ~ 3 x 10 18 eV (j25h. A 
single power-law injection source, modified by acceleration and transport 
effects, provides the simplest solutions. 

A comprehensive model of HECRs therefore seems achievable if a source type 
that injected power-law distributions of relativistic CRs were found in star- 
forming L* galaxies such as the Milky Way, with injection episodes frequent 
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and energetic enough to power the HECRs. For reasons reviewed in Refs. 
and (26), GRBs offer the most likely solution. 



2.1 Local GRB and Super-GZK Emissivities 



We first address the question of the local (z < 1) time- and space- averaged 
GRB luminosity density (or emissivity) £grb,xh (ergs Mpc -3 yr _1 ) and the 
local luminosity density Egzk required to power the super-GZK (> 10 20 eV) 
CRs. Pre-Beppo SAX estimates (|l||2j) concerning the X/7 emission from GRBs 
found 



with e 44 « few. This exceeds the local (z < 1) GZK emissivity eqzk — 
u/ttjyx ~ 6 x 10 43 w_2i ergs Mpc~ 3 yr^ 1 required to power > 10 20 eV CRs, 
where 10 _21 tt_2i ergs cm -3 is the observed energy density in > 10 20 eV CRs 
(u_2i = 0.5 for HiRes and w_2i — 2 for AG AS A), and the photo-pion energy- 
loss timescale for a 10 20 eV proton is 140 Mpc/c (|27l). A better estimate that 
corrects for the energy- dependence of implies Egzk — u/t^ ~ 10 44 m_ 2 i 
ergs Mpc~ 3 yr _1 , comparable with the GRB emissivity egrb,x/'y, so that the 
required luminosity density to power the super-GZK CRs is in coincidence 
with that available from GRBs. 

Based on the study by Bottcher and Dermer flljjh of BATSE statistics in the 
external shock model, a value of £44 = 4 was derived (0) using an analytic fit to 
the SFR derived from Hubble Deep Field measurements (|2al . This emissivity 
includes, in addition to the photon luminosity that is dominated by hard X-ray 
and 7-ray emission, an inefficiency factor for production of photons due to the 
radiative regime implied by the fits to the data 

(0), 

and a further correction 
for dirty and clean fireballs (|29l). which have since been discovered (|3oh. 

Removing the inefficiency factor from the estimate of Ref. (0) gives the local 
GRB emissivity in the form of nonthermal photon radiation to be £44 = 0.6. 
Vietri, de Marco, and Guetta f)3lh argue for a value of £44 = 1.1. Given the 
various uncertainties, we take the local GRB emissivity in the form of hard 
X-rays and 100 keV - MeV 7-rays from GRBs to be 



tGRB,xh ^ 10 44 e 4 4 ergs Mpc 3 yr 



1 



(1) 



£44 = 1 . 



(2) 
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2.2 Beaming and GRB Rate 



A second issue to be considered is the source rate of GRBs in our Galaxy. 
Beaming breaks in the optical light curves of GRBs indicate that the most 
apparently luminous GRBs are highly collimated. Frail et al. (Jiih argue that a 
typical GRB has a beaming factor of l/500 t/l of the full sky, so that GRBs are 
in actual fact 500 times more numerous, and l/500 th as energetic on average, 
as the measured energy releases imply. This means that most GRBs typically 
release ~ 5 x 10 50 ergs in X/7 emission. 

The local density of L* spiral galaxies, of which the Milky Way is representa- 
tive, can be derived from the Schechter luminosity function, and is ~ 1/(200- 
500 MpcJ}_ S H). The BATSE observations imply ~ 2 GRBs/day over the 
full sky (34). Due to beaming, this rate is increased by a factor of 500/soo, 



where /500 ~ 1- Given that the volume of the universe is ~ 47r(4000 Mpc) 3 /3, 
this implies a rate per L* galaxy of 

(200 - 500) Mpc 3 /L* 700 Cf __ 
PL * * f (4000 Mpc)3 ^ X 500/500 X SFR X Kft 

« (200 - 500) ^ x 0.33 x 500 /500 x (*™) x K FT 

L* Gpc d yr 1/8 

« (1 - 3) x 10- 4 (^) x (^) x f 50Q yr- 1 . (3) 

The factor SFR corrects for the star-formation activity at the present epoch 
[SFR(z = 0) = (l/8)SFR(z = 1)] (see Section 4), and the factor K FT ac- 
counts for dirty and clean fireball transients that are not detected as GRBs 
(jilt ). This estimate is in agreement with the result of Ref. (31) that excludes 



short GRBs. Thus a GRB occurs about once every 3-10 millennia throughout 
the Milky Way, or at about 10% of the rate of Type Ib/c SNe. 

We write Ecr = fcR^GRB,x/-y f° r the local luminosity density of CRs injected 
by GRBs, where fcR is the nonthermal baryon-loading factor by which the 
emissivity injected by GRBs in the form of hadronic CRs exceeds the emissiv- 
ity inferred from direct observations of the X/7 emission from GRBs. In the fol- 
lowing, we take £ 44 = 1 (eq. 2). From eqs. (3) and (1), we find that the average 
apparent isotropic energy release per GRB is = 10 53 (fcR + £44) /(K FT /3) ergs, 
and the smaller actual mean energy release per GRB due to beaming and the 
consequent larger number of GRB sources is = 2x 10 5O (/c*i?+£ ; 44)/[/5oo(-^'FT/3) 
ergs. 

For our canonical model GRB used to fit CR data near the knee, we take 
Ecr = 10 52 ergs, corresponding to }qr = 50. A value of fcR ^ 50 is implied 
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by the data fits if this model for HECRs is correct (Section 5). Hence, in a 
unified model for HECRs from GRBs, the GRB blast wave must be strongly 
baryon-loaded with f C R ^ 50, at least during the prompt phase of the GRB 
when the acceleration is most rapid. 



3 Galactic Cosmic Rays from GRBs 



Observations indicate that GRB sources are located in the disks and star- 
forming regions of galaxies undergoing active star formation, such as the Milky 
Way. The relativistic ejecta in the GRB explosions accelerate and inject CRs 
into the ISM of the GRB host galaxy in the form of a power law to the highest 
energies, though possibly with a low-energy cutoff to the accelerated proton 
spectrum. In relativistic blast waves, this cutoff is expected at ~ T 2 GeV 
energies, which could easily reach TeV - PeV energies for typical GRB blast 
waves with T ~ 10 2 - 300. UHECR acceleration and injection probably occurs 
during the 7-ray luminous phase of a GRB, which is on the order of minutes 
to hours. Acceleration and injection of lower energy CRs might operate on the 
Sedov time scale, which could exceed thousands of years. These acceleration 
times are still short compared to the times for particles with the corresponding 
energies to diffuse a distance comparable to the disk scale height, as can be 
shown from the diffusion properties of the Galaxy derived below. A GRB 
source can therefore be treated as an impulsive source of CR injection. 

In this section we present a simplified propagation model for CRs in the disk 
and halo of the Galaxy. 



3. 1 Diffusion Mean Free Path 



We assume that the Galaxy's disk and halo magnetic field consists of a large- 
scale field of mean strength B = B^ G fiG on which is superposed a spectrum 
of MHD turbulence. The diffusion mean-free-path A for CRs diffusing through 
this field is a function of the Larmor radius 

r L = J J~ sin a = — — = 3.12 x 10 12 — L- cm = — ^- pc , 4 
ZeB QB ZB^ ZB^ P ' y ' 



where 7 = 10 6 76 is the Lorentz factor of the relativistic CR proton or ion 
with atomic mass A, /3 — a/1 — 7 -2 , and the pitch angle a is set equal to tt/2. 
We assume isotropic turbulence, though it is straightforward to generalize the 
treatment for particle-scattering properties that differ in directions parallel 
and transverse to the plane of the disk. 
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We consider a very simplified approach to CR transport where particles diffuse 
via pitch-angle scattering with resonant MHD turbulence; see, e.g., Refs. (J35I 
36l : 37: liif). and Ref. f|39r ) for a detailed treatment. Let w(k)dk represent the 



differential energy density of isotropic MHD turbulence with wavenumbers 
between k and k + dk. First consider the case where the MHD turbulence 
spectrum is described in the inertial range by a single power-law function 
kw{k) = Wo{k/ki) l ~ q H(k; ki, k 2 ), where H(x;a,b) = 1 when a <x <b and 
H(x; a, b) = otherwise, and q is the spectral index of the wave spectrum. 
The index q = 5/3 for a Kolmogorov spectrum of turbulence, q = 3/2 for a 
Kraichnan spectrum of turbulence, and w tot = w /(q — 1) is the total energy 
density of MHD turbulence when q > 1 and ki 3> k\. Note that this description 
of the wave turbulence is an extreme oversimplification, insofar as we do not 
consider different wave modes and helicities, or distinguish between forward- 
and backward-propagating waves. 

The ratio of the MHD wave energy wtot to the energy density Ub = B 2 /87r 
of the large-scale component of the Galactic magnetic field is denoted by £i, 
so that kw(k) = (q - 1)^V B (k / 'k 1 ) 1 - q for k x < k < k 2 . The value of k x 
corresponds to the inverse of the largest size scale on which turbulence is 
injected. 

We make the ansatz that the diffusion coefficient of a CR with Larmor ra- 
dius tl is inversely proportional to the the energy density in gyroresonant 
waves with wavenumber k ~ r^ 1 (jUJ. We therefore have A = r^UB/kw(k) = 
kl~ q rl~ q /^{q - 1) = b^rl^/^iq - 1) oc 7 2 " 9 , giving 

a(pc) = r^r(^r) 2 - q = /2 , (5) 



ZB„ 



for q = 3/2 



where k\ = 1/bi = b 1 pc 1 . 

Equation (5) gives the mean-free path A of relativistic (7 ^> 1) CRs that diffuse 
due to gyroresonant interactions with large- wavenumber {k > kx) turbulence. 
When r^ 1 < k%, or > b pc pc, CRs gyroresonate with small-wavenumber 
(k < ki) turbulence. If the knee in the CR all-particle spectrum at — 3 
PeV is due to propagation effects, then b pc = 1 when B^ G = 3 for the dominant 
CR proton species. The break energies in the spectra of other CR ions are also 
given by the relation r L = fcf , so that the knee energies for the different CR 
ions occur at energies 

E z (PeV) = A l6 ZB^ G b pc . (6) 
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The fit to the data described below uses b pc = 1.6 and B^q = 3. If the value 
of £1, corresponding to the ratio of the energy density k ^ ki turbulence to 
Ub is a few percent, then A ~ ha for « 3 PeV CR protons, where /i^ is the 
scale height of molecular gas or massive stars in the disk of the Milky Way. A 
break in the all-particle spectrum due to a change in the propagation mode 
occurs at ~ 3 PeV, where A ~ ha, with a pattern of break energies Eb r oc Z 
for the ionic constituents of the CR all-particle spectrum. 

Although a CR proton with energy Ex has ~ pc <C h^, its mean-free- 
path \{E K ) ~ h d ^> r L (E K ) because pitch-angle scattering results only from 
gyroresonant interactions with MHD waves that constitute a small fraction 
of the total magnetic field energy density. CRs injected with E 3> Ek will 
diffuse through the Galaxy's disk and halo through gyroresonant pitch-angle 
scattering with MHD turbulence carried by waves with wavenumbers k <C k\ 
that gyroresonate with particles with E 3> Ex. This suggests that we consider 
a more general wave spectrum over a larger range of k given by the expression 

kw{k) = A [f H (^y-^H(k; k , h) + {^) l - q H{k; k h k 2 )\ , (7) 



where qn is the index of the small wavenumber MHD turbulence, and fu 
normalizes the energy densities of the small (k < ki) and large (k > ki) 
wavenumber turbulence spectra aXk = k\. The normalization condition ^Ub = 
J °° dkw(k) implies A = t^Us/g, from which the normalization coefficient g 
can be easily derived. 

The wave spectrum from eq. (7) is plotted in Fig. 1 by the solid lines. The pa- 
rameters are k = 1/100 pc -1 , k x = 1/1.6 pc" 1 , fu = 0.2, £ = 0.1, qn = 3/2, 
q = 5/3, and B = 3/iG. Sixty-five percent of the turbulence energy density 
is in the small wavenumber component, so that £i = 0.035. The sharp edges 
in this wave spectrum introduce unphysically sharp features in the spectra of 
particles for this simple model of transport. A more physically plausible tur- 
bulence spectrum that retains the features of eq. (7) but smooths the particle 
distribution function, as would be physically expected from 3-wave interac- 
tions and turbulence-energy diffusion and cascading (|4lh. is given by 

kw(k) = /^(±)i-^A A (£)[l - A A (A)] + A (^[l - A A (|)], ( 8 ) 

where the smoothing function A A (x) = exp(— x" A ), and we let k 2 — > oo. 
Eq.(8) reduces to eq. (7) in the limit A > 1. Fig. 1 shows a plot of eq. (8) 
with A = 1.5 and, for comparison, A = 6. The subsequent fits to the data use 
the MHD wave spectrum, eq. (8), with A = 1.5. 
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The ansatz A = r L [UB/kw(k)] with k — > r L : therefore implies 

A ^ — (Q) 

V ^(^O^^AaMOII - A A (^o)] + (Jlfci)'" 1 !! - A A (rLfci)] ' 1 J 



In Fig. 2, we plot A and r^ as a function of total energy E for CR protons 
and Fe nuclei. The turbulence spectrum is given by Fig. 1 with A = 1.5. Note 
that A = 10 kpc when the CR proton energy E 4 x 10 17 eV. This energy is 
near the onset of quasi- rectilinear propagation in a Galactic halo of size ~ 10 
kpc, and we expect a cutoff in the CR proton flux from a recent GRB in this 
energy range whenever we are not in the beam of a GRB jet (which is nearly 
always true because of the extreme improbability of such an occurrence). The 
energy-dependent cutoff will be oc Z, as observed. The actual situation is of 
course more complicated due to the magnetic field gradient from the disk to 
the halo and into metagalactic space. 



3.2 Three Dimensional Isotropic Diffusion 



We assume that a GRB injects a total energy Ecu = 10 52 E 52 ergs in the 
form of cosmic rays. For impulsive injection, the number injection spectrum 
of CRs with atomic charge Z and atomic mass A is assumed to be given by 

Nz,a{i) = c z ,aK((3^)~ p for ((3^) min < < {(3^) max , and we assume that p 
is the same for all (Z, A). When (f3^) min < 1, {f3^) max > 1, and 2 < p < 3, 

K = [1 + ^t]- 1 {p ' 2)Ec ' r . (10) 
L 2(3 -pY m p c 2 J2z,A Ac z,A V ; 



If 1 < (Pi)mm = Imin < (Pl)max, 2 < p < 3, 

(v- 2)E r Rl p ~ 2 

r-u \y ^J-^OH Imin (11^ 

mpC 2 J2z,a Ac z ,a ' 



Cosmic rays are assumed to diffuse isotropically in the interstellar medium, 
with the initial anisotropy from injection by the GRB jet quickly washed 
out. For energy- dependent diffusive propagation, the spectral number density 
of relativistic CR ions with charge Z and mass A measured at a distance 
r = 500r 500 pc away from an impulsive source of CRs from a GRB that 
occurred a time t earlier is 

nzAr,r,t) = ° Z 'tn2 exp[-(r/r di/ ) 2 ] . (12) 
n r dif 
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Here the diffusion radius r<&/ = 2^jD{pf)t = 2y Act/3, and eq. (12) assumes 
that CRs suffer no significant energy losses during transport (see Ref. (fiih for 
the more general case). 

Inspection of eq. (12) shows that when the observer is within the diffusion 
radius, that is, when rCr^j, then 

7 -p 7 -p-| (2-9) _ f T v ~^ , for q = 5/3 

Wz,a(7; t™, *) oc — oc -p oc i 3/2 x < . (13) 

^ I 7 - p -3 , for q = 3/2 



The measured spectrum from a burstlike source is therefore steepened by |(2 — 
q) units, where the diffusion coefficient D oc A oc 7 2-9 (0). (By comparison, 
the spectral index is steepened by (2 — q) units for continuous injection sources 
of CRs.) Eq. (13) shows that energy- dependent diffusion due to pitch-angle 
scattering with a Kolmogorov and Kraichnan spectrum of MHD turbulence 
steepens the injection spectrum of an impulsive source by 0.5 and 0.75 units, 
respectively. If the disk and halo magnetic fields of the Galaxy support a two- 
component MHD turbulence spectrum with indices q = 5/3 and qu = 3/2, 
then an injection spectrum with p = 2.2 will be steepened to a measured 
spectrum 7^,4(7; r, t) oc 7 _s , with s = 2.7 at E <C Ek and s = 2.95 at 
E ^> Ek- Because these indices are similar to the measured CR indices below 
and above the knee energy, we adopt this simplified model for CR transport. 
Speculations concerning the origin of such a turbulence spectrum are deferred 
to Section 7. 



Fig. 3 shows the CR proton flux measured at Earth at different times after 
a GRB explosion located 500 pc from Earth. The low-energy cutoff energy is 
given by (^j) m i n <C 1, so that most of the energy is deposited by ~ 1 GeV/nuc 
CRs. When {(3^) m i n = {(3"i) m i n ,\ 1, the CR intensity is increased at energies 
above the low-energy cutoff by the factor [(/3 , y)min,i] p ^ 2 ■ 

The flux of relativistic CRs is calculated from the expression 

10 4 c 

Iz,a(E)(™ 2 s sr GeV)- 1 = ^ E{Q(N) in z , A ( T ,r,t) , (14) 



where E = 'jAm p c 2 and riz t A is in cgs units. Propagation effects arising from 
the combined disk and halo MHD turbulence spectrum produce a break in the 
spectrum at E ~ 3 PeV. The index of the CR number fluxes above the knee 
is ~ 3.0, for the reasons given above, until propagation effects at the "second 
knee" at E pa 4 x 10 17 eV start to soften the spectrum so that the metagalactic 
component begins to make a dominant contribution. 
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The observer is within the diffusion radius when ^^(7) = 2JD(pj)t 3> r, or 
when ct ^ r 2 /X. If we are considering CRs with energies below the knee, then 
these particles scatter with the large wavenumber k > hi, q — 5/3 turbulence, 
so that A = gk{~ q r1~ 9 /£. We write ct = 3.07 x 10 5 tM yr pc, where the observer 
is irradiated by a CR flux from a GRB that took place tuyr Myr ago. At 
early times following the GRB, the low-energy portion of the measured CR 
spectrum from a GRB source is exponentially attenuated due to the slower 
diffusion of lower energy CRs. The observer will see the maximum flux from 
CRs with energies E = Am p c 2r y when 

t Myr - 0.6(^)r 5 2 00 -^(^)- 1/3 - O.OMr^-^)-^ . (15) 
9 bp'c Z± VG ZB^g 



This expression uses the parameters for the model fits discussed below, with 
£ = 0.1 and g = 4.3, thus showing that the maximum flux of CR protons 
near the knee energy reaches an observer 500 pc away m 14000 years after the 
GRB, in accord with Fig. 3. Eq. (15) shows that the low-energy cutoff jcutoff 
of the CR proton flux in Fig. 3 evolves according to the relation Jcutoff oc 
r 1/(1-3/2) £-1/(2-9) _ N ote ^at eq. (15) is only valid when t > r/c. 

Fits to the CR data near the knee are compatible with a 500 pc distant GRB 
releasing 10 52 ergs in cosmic rays if the GRB took place ~ 10 5 yrs ago. This 
implies that the 0.1 - 100 PeV CR fluxes were up to 100 times brighter one or 
two hundred thousand years ago than they are today, though no test of this 
implication suggests itself. 



3.3 Fits to the KASCADE Data 



We now apply this model to fit the preliminary KASCADE data announced 
in 2001 (H; El E]) . The points in Fig. 4 show KASCADE data and error bars for 
CR protons, He, Carbon, Fe, and the all-particle spectrum in panels (a) - (e), 
respectively. Fig. 4(d) also shows the reported la boundary to the estimated 
systematic errors for CR Fe. 

The fits to the KASCADE data shown in Fig. 4 use the propagation model 
previously described, with a GRB source of CRs at a distance r = 500 pc 
that exploded 210,000 yrs ago. The change in the diffusion and propagation 
properties leads to a break in the CR particle spectrum at break energies 
oc Z (eqs. [5] and [6]). Because of the preliminary nature of the data and the 
potentially significant systematic errors that could still remain in the early 
analyses, we did not perform a rigorous fit to the data, but instead adjusted 
the wavenumber k\ and the compositions of the different ionic species until a 
reasonable fit to the data was obtained. With B = 3/iG, the best value for the 
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Table 1 

Cosmic Ray Abundances Relative to Solar Photospheric Abundances 1)43 ) 



Ion cz,a Enhancement 





1.00 


1.0 


fHe 


0.098 


1.6 


12p 
6 U 


3.63 x 10~ 4 


50 


16q 
8 U 


8.51 x 10~ 4 


£ 5 


56 rv, 
26 re 


4.68 x 10~ 5 


20 



energy-dependent break was obtained with k\ = 1/1.6 pc in the spectrum of 
turbulence. 



The Anders- Grevesse Solar photospheric composition (|43l). and the (energy- 
independent) metallicity enhancements, compared to the Anders-Grevesse com- 
positions, that were used to fit the data, are listed in Table 1. The strong 
enhancements by a factor of 50 and 20 for C and Fe, respectively, may be 
possible from highly-enriched winds that could be peculiar to a GRB stellar 
progenitor, or if there was an earlier supernova as in the supranova model. 
Note that O cannot be strongly enhanced if we are to maintain a good fit to 
the all-particle spectrum. A more detailed discussion about implications from 
composition will depend on fits to final analysis of KASCADE data (44). 



As can be seen from Fig. 4, this model gives a reasonable fit to the data 
for the CR ion and all-particle spectra around the knee. Thus we argue that 
propagation effects from CRs injected and accelerated by a single GRB are 
responsible for the shape of the CR spectrum between « 1 and 100 PeV. 



This model must not overproduce CRs at energies below the knee. Fig. (5) 
compares CR observations with a model that fits the KASCADE data and has 
an injection momentum at ((3j) m i n = 0.01. This model weakly overproduces 
the CR Carbon ions, though it is not a serious discrepancy given the uncer- 
tainty in the KASCADE data which determines the CR Carbon composition. 
A strong overproduction would imply a low-energy cutoff in the energy of the 
injected cosmic rays, or an earlier GRB. A model with a low-energy cutoff at 
10 5 GeV would certainly be consistent with the medium-energy cosmic rays. 
In the proposed scenario, SN contributions to CR production make up the 
difference between CRs produced by a single GRB source and the measured 
CR fluxes below ^100 TeV/nuc. The transition from multiple SNs source to, 
primarily, a single source seems to occur at energies of ~ 1 - 100 TeV/nuc. 
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3.4 Probability of Nearby GRB and CR Anisotropy 



We estimate the probability that we find ourselves in the midst of a CR bubble 
formed by a recent, nearby GRB using the relation P ~ (r di f /r MW ) 2 NQ RB . 
The quantity (raif/rMw) 2 measures the relative area of the Milky Way's disk 
covered by a CR bubble with age t, and Nqrb is the number of GRBs that 
would have taken place during t. We normalize the rate of GRBs in the Milky 
Way to N M GRBs per millennium, and N M « 0.1-0.3 (Section 2). Thus 

NgRB = NcRBt = 10 3 iV M t M yr- 

Using model parameters, r^/ = 2 J Act/3 = AM\^ r (A^§ / Z B 1//6 kpc. Most 
of the active star formation in the Milky Way takes place within the spiral arms 
found within r MW = 15r 15 kpc. Thus P ~ 106(t Myr /r 15 ) 2 N M (A'j 6 /ZB flG ) 1/3 , 
and the probability is very unlikely when P(tMyr) <C 1, or when t Myr ^ i Myr = 
O.lrisA^ 1 ^ 2 (A^q/ZB^g)^ 1 ^ 6 . The model fits use a time of 0.21 Myr since the 
GRB exploded to fit the CR data, so there is a reasonable probability for this 
event to have taken place, even with A^m ~ 0.1. 

CRs from a GRB event that occurred much earlier than 0.2 Myr cannot fit the 
data unless the total energy release in CRs from a single GRB is increased, 
which makes severe demands on GRB models (see Fig. 3). For the estimated 
values of the diffusion coefficient, CRs from a GRB that occurred much later 
than « 1 Myr also makes severe energy demands on the GRB source. Thus 
a GRB event occurring a few hundreds of millennia ago represents the most 
probable situation. 

The anisotropy u = ilmax ~ I m in)/{I m ax + I mm) = (3D/cn ZA )(dn ZA /dr) = 
(^/nz,A)(dn>z,A/ dr) ()47l ). Substituting eq. (12) gives 

W = ^ = ^- - °-^%. (16) 

r dif 2 Ct tMyr 



The interesting point about eq. (16) is that u is independent of energy for this 
diffusion model. 

Analyses of the arrival directions of lar ge n umbers of CRs show anisotropics 
near 3 PeV at the w 0.15 ±0.05% level (£353). When E is between 100 GeV 
and 100 TeV, anisotropics of ~ 0.1% are found. Schlickieser (|39h discusses 
some sources of uncertainty in anisotropy measurements. An anisotropy of 
w 1% is implied from eq. (16) for tu y r = 0.21. If an anisotropy below ~ 0.2% 
is confirmed, then a number of implications follow. Either we are located 
near a rather recent GRB, which could be unlikely, or the CR energy release 
from GRBs is larger than given here. Moreover, there could be contributions 
from a second GRB that would help isotropize the flux. Indeed, CR contribu- 
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tions from a multitude of weaker SNe that do not host GRBs could also help 
isotropize the CR flux. We discuss the relative contributions of SNe to the CR 
spectrum in Section 7. 



4 Ultra-High Energy Cosmic Rays from GRBs 



In our model, the UHE component of the CR spectrum is produced by extra- 
galactic GRBs. Throughout our calculations we take the local GRB luminosity 
density in CRs to be e C r = fcR£GRB,x/~, where e GRB ,xh is S iven in ec l s - i 1 ) 
and (2), and Jcr is the nonthermal baryon- loading fraction required by a GRB 
model of HECRs. Our fits to the data imply values for fc R - 

We assume that the GRB cosmic rate-density evolution, in comoving coordi- 
nates, follows the SFR history (28) derived from the blue and UV luminosity 
density of distant galaxies, with an analytic fitting profile given by 

" w = " (0) (TT^T^TT^ (17) 



(|50f). where n{z) is the GRB comoving rate density at epoch z. To accom- 
modate present uncertainty in the true SFR evolution we take two models 
with a x = 0.005(0.0001), a 2 = 3.3(4.0), and a 3 = 3.0(3.0) in line with the 
extreme ranges of optical/UV measurements without (with) dust extinction 
corrections. The lower optical/UV curve is arguably a lower limit to the SFR 
evolution, while a much stronger evolution of the SFR is found by Blain et 
al. (|51l ) after correcting for dust extinction. To study the significance of this 
uncertainty we consider a "lower SFR" and "upper SFR" in each of our model 
fits. We plot these SFR histories in Fig. 6 in comparison with an oversimplified 
n{z) = n(0)(l + z) A evolution used by other authors (J59t lei} ), who perform 
similar UHECR calculations. 



Both our upper and lower SFR evolution differ significantly from the evolution 
oc (1 + z) 4 , when z 2. For > 10 20 eV CRs, this difference is irrelevent as all 
of these CRs must have orginated from z< 1 due to severe CR attenuation 
during propagation, but the differing SFRs have an important effect on the 
predicted CR flux at energies ^ 10 18 eV. 

Results from the Wilkinson Microwave Anisotropy Probe (|61[ ) and the High-z 
Supernova Search Team (Jfj^) favor a A-cold dark matter (ACDM) cosmology 
where Q tot = 1.0, Q m = 0.3, and Qa = 0.7, assuming a dark energy equation 
of state with w = —1 and h = 0.7. We use this set of cosmological param- 
eters throughout. Our treatment of the UHECR flux calculation, including 
attenuation during propagation, is presented in Appenix A. 
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4.1 KASCADE and HiRes Data Fits 



Our fits to the KASCADE and HiRes-I and HiRes-II Monocular data cover 
the energy range from « 2 x 10 16 eV to 3 x 10 20 eV. The low energy CRs are 
from galactic GRBs (Section 3) and the high energy flux is a superposition 
of extragalactic GRB sources. The region where the combined flux changes 
from halo-dominated to extragalactic-dominated is near E^°. We carry out 
8 different calculations where in each case our minimum-x 2 routine fits the 
halo and extragalactic normalizations and finds the best-fit value for E^°. 
The value of E^° in turn fixes parameters in the halo propagation model. 

We take characteristic values for the maximum cutoff energy from acceleration 
in relativistic shocks of E max = 10 20 eV and 10 21 eV. The spectral indices we 
consider are p = 2.0, where the CR energy is injected equally over all decades in 
energy-space, and p = 2.2 where optimal fits are found in our Galactic diffusion 
model. The upper and lower SFR evolution cases are studied in turn. For each 
case our model has three free parameters: 1) the CR halo-component cutoff 
energy E^°; 2) the local CR luminosity density of extragalactic GRBs; and 3) 
the relative intensity of the galactic halo CR component to the extragalactic 
component. We also consider two values of the minimum energy E m i n = 10 9 eV 
and 10 14 eV for CR injection, which yield identical fits but affect the overall CR 
energy requirements of a typical GRB. For a soft p = 2.2 spectrum, the baryon 
loading f CR is a factor « (10 14 /10 9 ) - 2 = 10 lower for the E min = 10 14 eV 
case than the E min = 10 9 eV case. In total, we fit 43 data points covering the 
highest 4 decades of CR energies. We estimate the quality of our fits with a 
reduced-minimum-^ 2 routine for 40 degrees of freedom. 

The transition region between the halo and extragalactic components falls in 
the range ~ 2 x 10 17 eV - 2 x 10 18 eV. Thus, our data set must extend to 
low-enough energies to make a determination of practicable. The lower 
energy data we fit is from the KASCADE all-particle spectrum (0), starting at 
10 16 35 eV and continuing to 10 1715 eV. The HiRes-II Monocular data set covers 
10 17 ' 25 eV to 10 19 - 60 eV, while the HiRes-I Monocular data covers the highest 
energies, from 10 18 55 eV to 10 20 10 eV. Our best-fitting model (marked with 
an asterisk in Table 2) was considered in the context of the full KASCADE 
data set down to energies of 10 15 eV. 

The interplay between our model and features in the data is apparent in the 
8 calculations shown in Figs. 7-9 and with the results listed in Table 2. The 
steepening in the data near the second knee (10 17,6 eV) is a feature which rules 
out cases with p = 2.0. This assumes that the intercalibration between the 
KASCADE and HiRes data is accurate, a conclusion that could be modified if 
reanalysis of the HECR data reveal calibration uncertainties. The extragalactic 
flux for p = 2.0 spectra fall well below the data at E C r < 10 18 eV leaving the 
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Table 2 

The results of various model fits to the KASADE, HiRes-I and HiRes-II Monocular 
data. The parameters p and E max give the injection index and maximum injection 
energy of CRs from extragalactic GRBs, respectively, for the upper and lower SFR 
curves given in Fig. 6. The high-energy cutoff of the halo component and the 

required baryon load fcR are derived by minimizing x 2 , for low-energy CR injection 
cutoff E m i n = 10 9 (10 14 ) eV. The best fit model is marked with an asterisk. 



Spectrum, p 
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halo component incapable of fitting the data both above and below the second 
knee. Our best fitting set of model parameters is with p = 2.2, E max = 10 20 eV, 
and the upper SFR history. This case puts the transition between galactic and 
extragalactic CRs in the vicinity of the second knee, consistent with evidence 
for a heavy-to-light composition change in this energy range (jo^h . The ankle 
energy (at ^ 3 x 10 18 eV) has a simple interpretation as a suppression of the 
UHECR flux from the photo-pair process (analogous to the GZK suppression 
from photo-pion production). 

The extragalactic contribution to the CR flux in the 10 17 ' 6 eV - 10 18,6 eV 
range is driven by the strength of the SFR evolution. The flux below 10 18 ' 6 eV 
is markedly higher in the upper SFR case because the UHECRs produced at 
redshift z are readily swept down to ~ 10 18 / (1 + z) eV by propagation effects 
(see Appendix A). A similar result was obtained (jiot E3) in the case when 
the cosmic evolution of GRB sources changes oc (1 + z) 4 , which approximates 
the upper SFR to z ~ 2, but then diverges from the upper bound to the SFR 
derived in Ref. (jolT ). hence predicting too high a flux at E C r < 10 18 eV. 

4.2 AG AS A Data Fits 

The reported AGASA flux above 10 20 eV provides a clue that UHECRs may 
violate the GZK-cutoff. We explore this possibility within the context of an as- 



17 



trophysical source with a hard (p < 2.0) spectrum and cutoff E max = 10 21 eV 
sufficient to produce a super-GZK flux. The results of five models are pre- 
sented in Figure 10 with fits to AGASA's highest nine energy bins. For hard 
CR-spectra GRB sources, the CR energetics are not a problem because the 
high-energy portion of the spectrum receives most of the nonthermal CR en- 
ergy. However, it is clear from our fits that if a super-GZK flux is verified by 
observations with the Auger Observatory, it calls for either a multi-component 
model of astrophysical sources (including nearby sources) or new physics in 
the form of top-down (for recent reviews, see Refs. (64; 65)) or hybrid (for 
some examples, see Refs. (jo^it l68h) scenarios. 



5 High-Energy Neutrinos from Gamma-Ray Bursts 

The prediction of our model that the energy output in relativistic protons 
and nuclei accelerated by relativistic shocks of GRBs can be much larger than 
the total power inferred from the observed X-ray/MeV radiation, f CR ^> 1, 
suggests a very important conclusion concerning the question of detectability 
of GRBs by forthcoming km-scale neutrino detectors like IceCube. 

A recent study (flih of this question shows that, at best, only from the brightest 
fraction of GRBs with neutrino fluence Q u ,tot i n the prompt phase exceeding 
10~ 4 erg cm -2 , would it be possible to detect multi-TeV neutrinos by IceCube. 
This is a general model-independent result which only assumes that the radi- 
ation fluence derived from the X-ray/MeV 7-ray measurements in the prompt 
phase of GRBs does not strongly underestimate the total radiation output of 
GRBs (that could be missed at higher energies). This result follows from the 
direct calculation of the number of muon neutrinos N v to be expected for a 
detector like IceCube, with an area Ad = 10 10 Ai cm 2 , assuming a spectral 
fluence of the neutrinos e$ v (e) = 10~ 4 $_4 erg cm -2 . For a spectrum of e^ u (e) 
with the spectral index ~ (i.e., with index p u ~ 2 for the differential num- 
ber fluence (j) u (e)), a simple derivation give N v ~ 1.2$_ 4 A 10 (1 + ^lne^ 1 ) for 
eu < 1, and N u ~ 1.2$_ 4 Ai / v /ei4 for ej4 > 1 (see (JUJ). Thus if the neutrino 
fluence is comparable to the X/7 photon fluence, detection of v ^ neutrinos 
could be expected only from very powerful GRBs, at a fluence level $_4 3> 1. 

Note that $_4 = 1 corresponds to the neutrino fluence ~ 2.3 x 10~ 4 erg cm -2 
integrated per each decade of energy. Note also the decline in the number of 
neutrinos to be expected per equal energy fluences per decade when e i4 > 1. 
This is explained by the change in the detection efficiency P Ufl of upward- 
going muon neutrinos at this energy (see (@)), which can be approximated 
as P utl = 10 _4 ei 4 , where % = 1 for €14 < 1, and x = 0.5 for eu> 1. Therefore, 
since a typical spectrum of neutrinos expected from GRBs due to photomeson 
interactions of HECRs extends well beyond 100 TeV, detection of even a single 
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neutrino at multi-TeV energies would become likely only if the total fluence 
in the would significantly exceed 10~ 4 erg cm -2 . 

Meanwhile, the total amount of energy released in the neutrinos (including u e ) 
is only about 40-50 % of the energy of secondaries resulting from photomeson 
interactions. Therefore at least the same amount of energy is released in the 
gamma-rays and electrons of multi-TeV energies produced in the secondary 
pion decays. For a blastwave Lorentz factor T ~ 100-300, the GRB radiation 
fields in the prompt phase are typically optically thick to photo-absorption 
for gamma-rays with energies above the GeV domain as is apparent from the 
curves in Fig. 11. Even if one assumes that the radiation in the prompt phase 
of a GRB is contributed only by several individual spikes, with a character- 
istic duration as large as t spk ~ 5-10 s, the source becomes transparent to 
gamma-rays below the GeV domain only in case of a collapsar GRB with 
Doppler-factor 5 300. Therefore most of the energy injected in multi-TeV 
gamma-rays will be efficiently converted to hard radiation in the X-ray to sub- 
GeV gamma-ray domains through the synchro-Compton photon-pair cascade 
developing in the relativistic shock/GRB source. 

Our knowledge of high-energy 100 MeV) emission from GRBs is limited to 
7 GRBs detected with the spark chamber on EGRET 

0) and the Milagrito 

detection of GRB 970417a (71). The average spectrum of four GRBs detected 
simultaneously with BATSE and the EGRET spark chamber from ^ 100 MeV 
to ~ 10 GeV is consistent with a —2 photon number spectrum, implying that 
in these cases the ^ 100 MeV fluence does not exceed the BATSE fluence by 
more than a factor of a few. However, the examples of the prompt and delayed 
high-energy emission from GRB 940217 ((72), the anomalous hard component 
detected with BATSE and the EGRET TASC from GRB 941017 Q, and 
GRB 970417a, which required ~ 10 times more energy in the TeV range than 
in the BATSE range, show that the high-energy behavior of GRBs is yet poorly 
measured and even more poorly understood. The baseline assumption made 
here is to require that the energy fluence in neutrinos not exceed the fluence 
observed in electromagnetic radiation at hard X-ray and soft 7-ray energies. 

Calculations of the neutrino fluxes expected from GRBs in case of baryon- 
loading factor f CR ~ 1, i.e., assuming equal energies for relativistic hadrons 
and directly accelerated electrons (which are assumed to produce the observed 
radiation), show that even in the case of very powerful bursts with radiation 
fluence at the level Q ra d ~ 10~ 4 ergs cm -2 , which happen only few/several 
times per year, the probability of detection of neutrinos by a km-scale detec- 
tor would remain hopelessly small in the framework of the collapsar model 
unless the Doppler factor 5 ^, 200 or the variability time scale t var < 1 s. 
A significant contribution to the target photon field density from a radiation 
component external to the GRB source/blob, as in the plerionic emission in the 



"supranova" scenario (|74j) , or from radiation scattered by the progenitor mate- 



19 



rial in the circumburst environment in the collapsar model, improves chances 
to detect neutrinos from GRBs (see Ref. flit ) for details). This is because for 
large Doppler factors and <& ra d ^ 10~ 4 ergs cm -2 , only a small fraction of the 
injected proton energy can be converted into secondaries through photomeson 
interactions of HECRs with internal photons. 

The current model for the origin of galactic high and ultra-high energy CRs 
from GRBs requires that the nonthermal baryon loading factor fen 1, 
and this prediction will be tested by the IceCube class neutrino detectors. In 
Fig. 12a we show the neutrino fluxes expected in the collapsar scenario from 
a burst at the radiation fluence level $ ra d = 3 x 10~ 4 ergcm -2 , calculated for 
3 values of the Doppler factor 5 assuming a GRB source at redshift z — 1 
(for h = 0.65). In order to demonstrate clearly the dependence of the neutrino 
fluxes on 5, here we fixed the baryon- loading factor at a value fen = 20 
for 3 values of 5. As in Fig. 11, for calculations in Fig. 12a we assume that 
the prompt emission is contributed by N sp k = 50 spikes with characteristic 
timescales t sp k — 1 s each, which define the characteristic size (in the proper 
frame) of individual spikes through R' spk ~ t spk 5/(l + z). 

The numbers of (muon) neutrinos to be expected by IceCube for 5 = 100, 200 
and 300 are N u = 1.32, 0.105 and 0.016, respectively. We should note, however, 
that for the assumed value of fcR, the calculated total fluence of neutrinos 
(both Ufj, and u e ) produced when 5 = 100 is & u ,tot = 7.2 x 10~ 4 erg cm -2 , i.e., 
by a factor 7.2/3 = 2.4 larger than the assumed radiation fluence. It means 
that in the light of the arguments given above, the maximum value of the 
baryon loading that could be allowed if the high-energy radiation fluence is 
less than the X/7 fluence for this particular case should be about 8-10, instead 
of 20. Consequently, the expected number of neutrinos for 5 = 100 should be 
reduced to ~ 0.6. On the other hand, in Fig. 12a the neutrino fluence for the 
case 5 = 200(300) is equal to <& v>tot = 1.4 x 10~ 4 (3 x 10" 5 ) erg cm" 2 , so this 
accommodates a baryon-loading increased from 20 up to fen — 45(200), with 
the expected number of neutrinos observed by IceCube N UtCorr ~ 0.23(0.16). 
If the radiation fluence at MeV - GeV energies is allowed to exceed the X/7 
fluence by an order of magnitude, a possibility that GLAST will resolve, then 
the expected number of detected neutrinos could be increased correspondingly. 

In order to understand the degree of stability of the predicted numbers of 
neutrinos against the model assumptions, and most importantly with re- 
spect to the assumed number of spikes which determine the characteristic 
size/compactness of individual "contributing blobs," in Fig. 12b we present 
the results of calculations of the fluxes assuming N spk = 10, correspond- 
ing to t sp k — 5 s. In this case we chose a higher Jcr = 40, anticipating a 
lower than previously considered fraction of energy extraction from the rel- 
ativists protons. The heavy solid and dashed curves show in Fig. 12b show 
the fluxes calculated for 5 = 100 and 300, respectively, in the collapsar sce- 



20 



nario. The corresponding total neutrino fluences are equal to 6.9 x 10~ 4 and 
1.3 x 10~ 5 erg cm -2 , and the relevant numbers of to be expected for Ice- 
Cube are N u = 0.94 and 6.4 x 10~ 3 . Correcting these numbers for the max- 
imum acceptable level of neutrino fluence assuming $ Vtto t ^ &rad, results in 
N^corr ~ 0.4 and 0.15 for 5 = 100 and 300, respectively, which are indeed 
close to the numbers derived above in Fig. 12a. Note that for the latter case, 
the implied baryon- loading factor would be as high as fcR — 900, which could 
still be marginally acceptable, as discussed in Section 5.1 above. 

The two thin lines in Fig. 12b are calculated in the framework of a "supranova" 
scenario for the same values of 5 = 100 and 300, assuming an external radiation 
field density in the "plerionic nebula" with age t s d = 0.3 yr (see Ref. (J19J) for 
details of calculations). The number of neutrinos to be expected when 5 = 100 
would be practically the same as for the collapsar model, with N VtCarr ~ 0.5. It 
is noteworthy, however, that for 5 = 300, the predicted probability of neutrino 
detection would be greatly increased relative to the collapsar model prediction, 
namely N U)Corr ~ 0.24, which would at the same time demand a significantly 
smaller value for the baryon-loading factor, with fen — 15 only. 

Although these numbers are still smaller than 1, these calculations leave a 
non- negligible probability for detection of 1 - 2, and hopefully up to a few, 
neutrinos from some of most powerful GRBs during a reasonable observation 
time even in the case of a "collapsar" GRB with 5 = 300. Note that these 
neutrinos would be mostly at energies > 100 TeV where almost no background 
neutrinos are to be expected from a given direction to any single GRB (19). 
Such a detection would be a confirmation of high baryon-loading in GRB blast 
waves, and would also suggest a very significant contribution of the accelerated 
hadrons in the observed hard radiation through secondaries produced in the 
photomeson interactions. 

Another prediction of a hadronically dominated GRB model is the possible 
detection of gamma-rays in the multi-GeV energy domain by GLAST at a 
stage near the prompt phase when the density of radiation fields would still 
be sufficient for extraction of a small fraction of the energy of accelerated 
protons, but when the secondary gamma-rays would not be suppressed by 
the photo-absorption process. In case of relatively close GRBs, with z <C 1, 
these gamma-rays could then be detectable also by ground-based gamma-ray 
detectors like VERITAS and HESS with energy thresholds ~ 100 GeV. 



6 Discussion 

We have proposed a model where HECRs originate from GRBs. The CR flux 
near the knee is assumed to result from CRs produced by a single GRB which 
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has occurred relatively recently and not very far from us in the Milky Way. 
These CRs propagate diffusively in the Galactic disk and halo. The simple 
diffusive propagation model developed in Section 3 implies that the measured 
CR flux results from the modification of the injection spectrum of an impulsive 
CR source due to transport through a magnetic field with a given turbulence 
spectrum. Using a turbulence spectrum harder at smaller wavenumbers and 
steeper at larger wavenumebrs, we have fit the 200 f KASCADE data for CR 
ion spectra between w 1 and 100 PeV, and explained the change in the all- 
particle spectra from p = 2.7 to p = 3.0. A GRB releasing « 10 52 ergs in 
HECRs, located ~ 500 pc away, and occurring k2x 10 5 years ago, provides 
reasonable fits to the KASCADE data. 

Our model of a single GRB source making CRs at energies through the knee 
of the CR spectrum bears some similarity to the single-source model proposed 
by Erlykin and Wolfendale (J751 ) to fit data near the knee. These authors argue, 
however, that propagation ( "Galactic modulation" ) effects cannot explain the 
constant rigidity break of the CR ionic species, whereas we employ a propa- 
gation model to produce that break. In this respect our propagation model 
treats rigidity-dependent transport as in the model of Swordv |45T) and builds 
upon the detailed study of Atoyan, Aharonian, and Volk (|42r ) for the spec- 
tral modification effects due to energy- dependent diffusive propagation of CRs 
from a single source. Our model explains, moreover, CR data not only through 
the first and second knees but also at the highest energies. 

The turbulence spectrum that fits the CR spectrum near the knee employs 
a Kraichnan spectrum at small wavenumbers and a Kolmogorov spectrum 
at large wavenumbers. Turbulence is thought to be generated at the small- 
est wavenumbers or largest size scales, with subsequent energy cascading to 
smaller size scales (76). It is interesting to note that there are two crucial 
length scales in our turbulence spectrum, namely k^ 1 ~ 100 pc, and fcjf 1 ~ 1 
pc. The generation of turbulence in the disk and halo of the galaxy at the 
larger size scale could be associated with halo-disk interactions (e.g., through 
the interactions of high- velocity clouds with the Galactic disk), which would 
deposit turbulence throughout the disk and the halo of the Galaxy on a size 
scale h d rs 100 pc. The smaller length scale is typical of the Sedov length scale 
for supernova explosions in the disk of the Galaxy. Indeed, SNe would generate 
a large amount of turbulence energy which could make a distinct contribution 
to the turbulence spectrum in this wavenumber range. 

The origin of the different indices of the two components of the MHD tur- 
bulence spectra at small and large wavenumbers could be related to the time 
available for the turbulence energy injected at the different size scales to cas- 
cade to larger wavenumbers. Medium-energy CRs with energies between 10 9 - 
10 14 eV/nuc will diffuse by gyroresonant pitch-angle scattering in response to 
MHD waves with k 3> 10/pc. The model turbulence spectrum at large values 
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of k is given by a Kolmogorov spectrum with index q = 5/3, as seen in Fig. 
(1). Because medium-energy CRs are thought to arise from a superposition 
of many SNe, we can treat their transport in the framework of continuous 
injection. As noted previously, the measured index of CRs from continuous 
sources is steepened by a factor 2 — q = 1/3. If medium-energy CRs, whose 
measured number intensity index is ~ 2.7 |77l). result from many SNe that 
produce CRs which diffuse through pitch-angle scattering in a spectrum of 
MHD turbulence with index q = 5/3, then it is necessary that the injection 
indices of these medium-energy CRs lie between w 2.3 and 2.4. 

This is a surprising result, because it is generally thought that the strong 
shocks in SNe accelerate and inject CRs with an injection closer to 2.0 than 2.4 
ffTih . Such a soft injection spectrum could be avoided if medium-energy CRs 
diffuse in a turbulence spectrum with index q = 3/2. In this case, the measured 
index is steepened by 0.5 units compared to the injection index. If the large-fc 
component with the steeper index is superposed on the extrapolation of the 
small-fc component to large wavenmbers, then at sufficiently large values of 
k, the turbulence spectrum will change from a Kolmogorov to a Kraichnan 
spectrum, as illustrated by the long-dashed line in Fig. 1. In this case, a 
continuous injection scenario implies that an injection index of 2.2 yields an 
observed CR spectrum with a 2.7 number index. A p = 2.2 injection index is 
in accord with expectations from nonrelativistic shock acceleration. Note that 
wavenumbers of 10 3 pc -1 are gyroresonant with ~ 30 TeV CR protons. 

An important feature of recent GRB studies is their association with SNe. 
CR acceleration at SN or GRB shocks is crucial for the production of CRs 
from the lowest energies, GeV/nuc, to the highest energies ^ 10 20 eV. The 
different speeds of the SN shocks in the different types of SNe ejecta ranging 
from relatively slow Type II (~ 3000 - 10000 km s _1 ) to marginally relativistic 
Type Ib/c is important to produce the full cosmic ray spectrum ( 
therefore predict that HESS and VERITAS, with their improved sensitivity 
and imaging, will detect 7-ray emission from supernova remnants at a low 
level unless that SN has also hosted a GRB. 

At energies E ^ E^° rs few x 10 17 eV, CRs stream out of galaxies to form the 
metagalactic CR component. We assume that GRBs evolve with cosmic epoch 
according to the SFR history of the universe, so that most of the UHECRs are 
produced at redshift z ^ 1 when the SFR is greatest. The UHECRs have their 
spectrum attenuated by photo-pion, photo-pair, as well as adiabatic energy 
losses which become important for particle energies below 2 x 10 18 /(l + z) eV. 
Our best fit was found to have spectral index p = 2.2 and E max = 10 20 eV. 
For these parameters, a slightly better fit was found for the stronger GRB 
redshift evolution ("upper SFR"). Stronger GRB evolution contributes more 
of the extragalactic CR flux over the range 2 x 10 17 eV to 2 x 10 18 eV than the 
lower SFR case, giving a best-fit value E%%£ ~ 2 x 10 17 eV. This is an effect 



HQ. We 
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of CR attenuation from photopion and photopair production processes where 
all of the CR flux produced at E ^ 2 x 10 18 eV (and 1) are swept to lower 
energies. Crucial for deriving these conclusions are the relative calibrations of 
the KASCADE and High-Res experiments. 

Our best fit to the UHECR data gives a measure of the local GRB luminosity 
density eqr required in CRs. We find fcR ~ 70 for E min = 10 14 eV and 
5 = 300. This implies that GRBs must be baryon-loaded by a factor ^ 50 if 
this model for HECRs is correct. The precise value of fen varies with Doppler 
factor according to fcR ~ 70(300/<5) a4 , for a p = 2.2 injection spectrum. 
Nonthermal baryon-loading factors fcR 3> 1 implies that GRBs should also 
be luminous in high-energy neutrinos. In section 6 we calculate the fluence in 
neutrinos and predict that ~ 0.1 - few neutrinos of energy 100 TeV-100 PeV 
may be observable in IceCube a few times per year from individual GRB 
explosions, which depends on the radiation fluence measured from GRBs at 
<; 100 MeV - GeV energies, which is not yet well known. 

In a collapsar model calculation |l9h. the predicted number of neutrinos de- 
pends sensitively on 5, because the density of internal synchrotron photons 
is much larger for smaller values of 5. The X-ray flashes (j3ol ) may be "dirty 
fireball" GRBs, which are similar to classical long-duration GRBs, though 
with larger (thermal) baryon loading and smaller average values of S (29). 
This suggests that the X-ray flashes may be more promising candidates for 
neutrino detection than GRBs with peak photon energies of the vF v flux at 
several hundred keV energies. 

The low statistics of current data at the highest energies make unclear the 
presence or absence of a super-GZK UHECR component. The HiRes data 
appears to be consistent with the GZK effect, thereby favoring an origin of 
UHECRs in "conventional" astrophysical sources. If the UHECR flux is found 
to exceed predictions of a GZK cutoff, then this could be evidence of new 
physics, could indicate the existence of super-GZK sources of cosmic rays 
within ~ 50 Mpc from Earth, or both. Our current interpretation of the data 
is made uncertain by results from the AGASA experiment (|20n which reveals 
a super-GZK flux, in conflict with the HiRes and Fly's Eye data (0) (see, 



however, Ref. (|24f). who argue that the discrepancy is not serious given the 
uncertainties in calibration and the statistical variance of CR fluxes measured 
at such high energies). 

The possibility of new physics associated with a super-GZK CR flux has gener- 
ated much interest. Scenarios for a super-GZK flux include the Z-burst scenario 



66f l. magnetic monopole primaries (|67t l68f ). and the decay of supermassive 



relic particles (|64j). (For a recent review, see Ref. (65).) Another approach to 
generate super-GZK fluxes is from astrophysical sources distributed within 
a GZK-distance from Earth, ^ 140 Mpc. The Auger detector should resolve 
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this question by measuring m 30 events per year above 10 20 eV. 

In conclusion, we have investigated the hypothesis that GRBs are the sources 
of HECRs. Our model provides a unified source-type to fit all the CR spectri, 
from some minimum CR energy produced in GRBs, which could be somewhere 
& 10 14 eV, to a maximum energy ^ 10 20 eV. The total CR flux that we 
calculate is a superposition of CRs originating from a past GRB in our Galaxy 
in whose HECR halo we inhabit (with a gradual transition at lower energies to 
CRs contributed from ordinary SN events with sub-relativistic blast waves), 
and ^ 3 x 10 17 eV CRs and UHECRs orginating in extragalactic GRBs and 
GRBs at cosmological distances. The HECR all-particle spectrum and our 
model fit to this data are shown in Fig. 13. 

A GRB model for HECR production requires that GRB blast waves are 
hadronically dominated by nearly 2 orders of magnitude. This could be related 
to the injection process in relativistic shocks and the large proton/electron 
mass ratio. This has important consequences for 100 MeV - GeV - TeV 7-ray 
emission and high-energy (^,100 TeV) neutrino detection. Allowed hadronic 
production in GRB blast waves will make 7 rays of comparable fluence as 
the neutrino fluence. Insisting that the high-energy radiation fluence from a 
photomeson cascade can only reach the level of the X/7 emission measured 
with BATSE (unless it is missed by detectors operating at larger energies), 
we have shown that J£ 1 neutrino could be detected coincident with a GRB 
at the fluence level ^3x 10~ 4 ergs cur 1 by IceCube. 

If GLAST shows that GRBs are much more fluent at 100 MeV energies 
than at X/7 energies, then in these GRBs one could expect a few high-energy 
neutrinos. Detection of even 1 or 2 neutrinos from GRBs with IceCube or 
a northern hemisphere neutrino telescope would unambiguously demonstrate 
the high nonthermal baryon load in GRBs. High-energy neutrino detection, 
especially from GRBs with bright 100 MeV - GeV emission components, or 
from GRBs with small Lorentz factors or large baryon-loading such as the 
X-ray bright GRBs, would provide compelling support for this scenario for 
the origin of the cosmic rays. 



We thank the referee for a helpful report. The work of S.D.W. was performed 
while he held a National Research Council Research Associateship Award at 
the Naval Research Laboratory (Washington, D.C). The work of C.D.D. is 
supported by the Office of Naval Research and NASA GLAST science inves- 
tigation grant DPR # S-15634-Y. A. A. acknowledges support and hospitality 
during visits to the High Energy Space Environment Branch of the Naval 
Research Laboratory. 



25 



A UHECR Attenuation and Flux Calculation 



The Larmor radius of a UHE proton is ~ 100 E^/B^q kpc for a proton of 
energy E = 10 20 E 2 o eV propagating through a magnetic field of strength 
B = B^ G 11Q. We approximate the UHECRs as propagating rectilinearly 
where their spectrum is modified by energy losses. Here we describe the 
three attenuation processes included in our propagation code: 1) red-shifting 
of the proton momentum, 2) pion production from proton-CMBR scattering 
(p 7 ~~ > 7r+ n or 71-0 P)i an d 3) electron-pair production from proton scattering 
on the CMBR (p 7 -> p e+ e~). 

The expansion of the universe causes a red-shifting of the proton momentum 

(A.l) 



1 dEJz) 1 dz 



E dt (1 + z) dt 

where for a ACDM cosmology 

dz 



df {l+z)H ^tt m (l + zf + tt A . (A.2) 



Energy loss from proton scattering on CMBR is given by 

1 dE v Jz) f , 71^(6,2) /" , / j- / \ . . 

_ E tft = J 27 2 e 2 y ^ 0-p 7 ( e »-)^P7( e r-) e r (A.3) 



(|80fl. where n 7 (e, z) is the black-body photon distribution of the CMBR with 
red-shift dependent temperature T(z) = 2.72K (1 + z), a P7 (e r ) is the cross- 
section for the process, and K p ^(e r ) is the inelasticity. The mean- free-path 
between scatterings is found from eq. (A.3) by setting K pi (e r ) = 1 for all e r . 



We approximate photo-pion production following Ref. (|81fl. modified to in- 
clude three step cross-section and inelasticity functions 

' 40 mb , for 150 MeV < e r < 220 MeV , 

(j p7 (e r ) = J 300 mb , for 220 MeV < e r < 450 MeV , (A.4) 

, 110 mb , for 450 MeV < e r , 
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and 



f 0.2 , for 150 MeV < e r < 220 MeV , 
0.2 , for 220 MeV < e r < 450 MeV , 
I 0.5 , for 450 MeV < e r . 



(A.5) 



With increasing e r , the three regions correspond to 1) direct pion production 
at threshold, 2) the enhanced cross-section at the A-resonance, and 3) multi- 
pion production in the diffractive scattering regime where photon coupling to 
the vector mesons po and uj dominate. In comparison with the attenuation 
calculations of Miicke et al. (jH) and Stanev et al. (27), we agree to ~ 2 (4)% 
at E p = 10 20 (10 21 ) eV. 

Our treatment of the photo-pair production losses follows that of Chodorowski 
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et al. (J83J). The greatest photopair energy- loss rate is found at E p ~ 3 x 10 
eV when z = 0, and pair-production dominates other energy losses over the 



range 2 x 10 18 eV £ E p £ 6 x 10 19 eV. 



Our UHECR flux calculation follows the formalism of Berezinsky and Grigor'eva 
(Jib! ) for diffuse sources. The flux from a volume element dV(z) = R z (z)r 2 drdVt 
observed at Earth is 



dJ(E)dE 



n{z)dV(z)F(E*,E max )dE* 
(1 + z)4nR%r 2 



(A.6) 



where n(z) is the comoving source density described in eq. (17) and Fig. 6, 
F(E, E max ) is the number of CRs produced per unit energy per unit time for 
an injection spectrum with power-law index p and exponential cutoff energy 
i? max , the Hubble length R = 4286 Mpc, and E* is the proton energy at 
generation which is subsequently measured with an energy E < E*. We use 
Ro = R(z)(l + z) and R(z)dr = cdt to obtain 

1 Zm f ax dt /E*\ 

J(E) = — F(E,E max ) J dz — n(z) (-) 

o 



dE*(E,z) 



dE 



(A.7) 



where the term \dE*(E, z)/dE\ transforms the energy interval from propaga- 
tion losses and we take z max = 4.0. Results of calculating eq. (A.7) for several 
different spectra and evolution models are found in Figs. 7, 8, 9, and 10. 
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Fig. 1. Wave turbulence spectrum used to model CR propagation in the Galaxy. 
The mean magnetic field is assumed to be 3 ^G. An idealized model is shown by 
the solid lines and, after smoothing, by the dotted and short-dashed curves for the 
smoothing parameter A = 1.5 and 6, respectively. The parameters for the model 
used in the fits are £ = 0.1, fjj = 0.2, ko = 1/100 pc _1 , k\ = 1/1.6 pc _1 , and 
A = 1.5. The extension of the small-wavenumber (k < k\) turbulence spectrum 
intersects the large- wavenumber (k > k\) spectrum at k ~ 1000 pc _1 , as shown by 
the long-dashed line. 
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Fig. 2. Larmor radius tl and mean-free-path A of CR protons (solid curves) and 
Fe nuclei (Z = 26, A = 56; dot-dashed curves) with total energy E in a magnetic 
field with mean strength of 3 /j,G. The wave turbulence spectrum given by Fig. 1 
with A = 1.5 is used to calculate A. 
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Fig. 3. Temporal evolution of the fluxes of CR protons (heavy curves) measured by 
an observer at Earth from a GRB that took place 500 pc away at earlier times given 
by the labels on the curves. The GRB was assumed to inject 10 52 ergs total energy 
in cosmic rays with an ionic composition given by Table 1 and a low-energy cutoff 
at (P^) m in = 0.01 (~ 50 keV). Shown by the thin solid curve is the all-particle flux 
at time t = 0.21 Myr used to fit the KASCADE data. Note that the fluxes have 
been multiplied by E 2 - 5 , where E is the CR particle kinetic energy. 
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Fig. 4. Data points show preliminary KASCADE measurements of the CR proton 
(panel a), He (panel b), Carbon (panel c), Fe (panel d), and the all-particle spectrum 
(panel e), along with model fits (dotted curves) to the CR ionic fluxes. In the model, 
a GRB that occurred 2.1 x 10 5 years ago and at a distance of 500 pc injects 10 52 
ergs in CRs, with ionic compositions given in Table 1. The CRs isotropically diffuse 
via pitch-angle scattering with an energy-dependent mean-free-path A in an MHD 
turbulence field given by Fig. 1. 
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Fig. 5. Contributions to the CR flux in the GeV/nuc - TeV/nuc range for a single 
source with a low-energy cutoff determined by the minimum injection momentum 
(l3^f)min = 0.01. Data is taken from the papers by Swordy et al. The excess 

flux from the model, especially for CR Carbon, suggest a low-enery cutoff in the 
CR injection spectrum from GRBs. 
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Fig. 6. The evolution of the rate density of GRBs as a function of 1 + z, shown in 
comparison to measurements of the star formation rate (SFR) inferred from mea- 
surements of restframe optical/UV emission from galaxies. The data, normalized 
with respect to a SFR of 0.09 M Q Mpc~ 3 yr -1 , are listed in order of increasing 
redshift and denoted for different measurements by the following symbols: filled 
hexagonsjE^) , open hexago ns JH^h . diagonal crosses (54), filled squares (jHH), open 
squares (|56l ) , and open stars ([571 ) . A SFR modification for dust extinction is given by 
the open triangles (|58l). We use the function given in eq. (17) with a = 0.005(0.0001), 
b = 3.3(4.0), and c = 3.0(3.0) for each of the models tested in our UHECR calcula- 
tions. The dotted curve, based upon optical/UV observations, is arguably a lower 
limit to the SFR evolution. An extreme correction for dust extinction (51) gives the 
stronger evolution shown by the middle, solid curve. For comparison, the dashed 
curve displays the relation n(z) = n(0)(l + z) 4 used by other authors 
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Fig. 7. Fits to UHECR data for model with injection spectral index p = 2.0. The 
fits to the K AS CADE (crosses), HiRes-I Monocular (squares), HiRes-II Monocular 
(stars) data assume various spectral cutoffs of the source, and lower and upper SFR 
histories and E max as labeled on the figures. Although the AGASA data (triangles) 
are shown, they are not included in the fits. The cutoff energies for the halo compo- 
nent E^°, reduced Xr> an d requisite local CR luminosity densities fcR are given 
on each figure. 
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Fig. 8. Upper-left figure gives a comparison of the upper versus lower SFR evolution. 
Attenuation of the CR flux originating at energies > 2 x 10 18 eV and z > 1 are 
swept by photo-pion and photo-pair processes down to lower energies. The increased 
number of CRs injected for the upper SFR yields a significantly larger flux of CRs 
in the range 2 x 10 17 eV to 2 x 10 18 eV. The remaining three figures have the same 
parameters as in Fig. 7, except that the injection index is now p = 2.2. 
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Fig. 9. Best fit to the Kascade (crosses), HiRes-I Monocular (squares), HiRes-II 
Monocular (stars) data assuming a spectral cutoff at the source of E max = 10 20 eV 
and using the upper limit to the SFR evolution. We also show the AGASA data (tri- 
angles) but do not include these in our fits. The cutoff energy for the halo component 
is E^a° = 10 17 ' 07 eV giving a Xr = 1-03 The requisite baryon loading factor for this 
fit is f C R = 746(70.3) for a low energy cutoff at the source of E min = 10 9 (10 14 ) eV. 
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Fig. 10. Best fits to the highest nine energy bins in the AGASA data (tiangles) 
3 x 10 19 eV - 3 x 10 20 eV for various hard spectra, p = 1.0, 1.2, 1.4, 1.6, 1.8. The 
resultant Xr an d requisite baryon loading fen for each case are shown on the figure. 
The KASACADE and HiRes data, although shown, are not included in the fits. 
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Fig. 11. Opacity r 77 for gamma-ray absorption in the radiation field inside in- 
dividual spikes of a GRB at z = 1, with the X-ray/MeV gamma ray fluence 
&rad = 3 x 10 erg cm -2 (with spectral flux peaked at 1 MeV) in the prompt 
phase with a duration tcRB = 100s, calculated for a collapsar GRB (a) (left): 
for 3 different Doppler factors 5 = 100 (solid curve), 200 (dashed curve), and 300 
(dot-dashed curve) , assuming that the prompt emission is contributed by = 50 
individual spikes with the mean duration t sp k = Is each; (b) (right): for 5 = 100 
(dot-dashed curve) and 300 (dashed curve), but assuming N sp f, = 10 spikes with 
the mean duration t sp k = 5 s. For comparison, here we also show contributions to 
r 77 due to the external radiation field in the case of "supranova" GRB with the as- 
sumed delay between the supernova and GRB events t s d = 0.3 yr; the 3-dot-dashed 
curve corresponds to the opacity due to the external photon field, and the solid 
curve shows the total the spike in the case of 5 = 100. The dotted curve 



shows the opacity across the entire "plerionic nebula" with a radius 
(see Ref. 1)1 9h for details of calculations in the "supranova" scenario). 



4.610~ d pc 
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Fig. 12. The fluences of muon neutrinos calculated for a collapsar GRB assuming 
(a) (top): 3 different Doppler factors 5 = 100, 200 and 300, the same GRB param- 
eters as in Fig. 11a, and a nonthermal baryon-loading factor fcR = 20; and (b) 
(bottom): Doppler factors 5 = 100 and 300, but for the case of t sp k = 5 s; for com- 
parison, here we also show by thin curves the fluences expected in the "supranova" 
model with the same parameters as in Fig. lib. 
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Fig. 13. The all-particle CR data from KASCADE, High-Res I and II, in comparison 
with the model result for the all-particle spectrum (dotted curve) from galactic 
(lower energy solid curve) and extragalactic (higher energy solid curve) GRB sources 
of high-energy cosmic rays. 
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